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Qv I In this work we derive certain topological theories of transverse vector 



c 
^ 



fields whose amplitudes reproduce topological invariants involving the inter- 
actions among the trajectories of three and four random walks. This result 
is applied to the construction of a field theoretical model which describes the 
K^ ^ statistical mechanics of an arbitrary number of topologically linked polymers 

Q^ ' in the context of the analytical approach of Edwards. With respect to previous 

o' 

^O I attempts, our approach is very general, as it can treat a system involving an 

Qv ' arbitrary number of polymers and the topological states are not only specified 

by the Gauss linking number, but also by higher order topological invariants. 

I. INTRODUCTION 

Topologically linked random chains are studied in connection with physical systems in 
which the topological constraints of one- dimensional objects play an essential role [|l|, [Q. 
Examples are vortex rings in fluids and dislocation lines in solids. To be concrete, the theory 
of random chains is applied here to the description of polymers in a good solvent 0- 0, 
where the random chains are subjected to excluded volume interactions and topological 
interactions |[T|, p- |0. The former take into account the effective repulsions experienced 
by the polymers, while the latter arise due to presence of stable topological constraints in 
the system. 

In a recent publication [^], the so-called analytical approach of Edwards ||^ has been 
extended to the case of N entangling polymers, mapping the statistical mechanical problem 



Oh! 



X 



of computing their free energy to a field theoretical problem. The model obtained in this 
way is a 0{n) field theory coupled to Chern-Simons (C-S) terms [|10| in the limit n -^ 0. At 
one loop approximation, one finds that the topological interactions tend to counterbalance 
the repulsive effects of the excluded volume forces. This effect has been indeed observed in 
nature, for instance in the DNA of some bacteria, which forms topologically entangled rings 
TT| . Moreover, one can compute the second topological momentum of two polymers exactly. 



The analytical approach of Edwards is however limited by the fact that the topological 
states of the system are distinguished using the Gauss linking number. The latter is a 
relatively poor topological invariant and describes topological interactions in which only the 



trajectories of two polymers are involved. Despite many efforts [|T2|^ IMIl) ^^^ inclusion in the 
treatment of topological random walks of more sophisticated link invariants, like for instance 
the Alexander and Jones polynomials [Q, |2^, has not yet been achieved. The reason is 



that higher order link invariants have no immediate relation to the physical conformations 
of the polymers in the space |]T^. A remedy invoked by many authors is the introduction 



of C-S field theories coupled to the polymer trajectories IQ, 0- Unfortunately, this is not a 
simple task. First of all, the amplitudes of non-Abelian C-S field theories contain topological 
invariants P3], but it is not clear how to use them to impose constraints on the configurations 



of the system. On the other side, it is not easy to find a regularization or a mechanism that 
suitably removes the spurious non-topological contributions arising in C-S field theories. For 
example, the introduction of a framing [^ complicates the integrations over the polymer 



trajectories to the extent that the mapping of the statistical mechanics of the system to a 
field theory is no longer possible [^ . 



To solve the above difficulties, we propose here interacting Abelian field theories contain- 
ing transverse vector fields. Such theories are topological but not gauge invariant. Of course, 
since the radiative corrections are not protected by a gauge principle, counter-terms may 
arise which are of a non-topological nature. To exclude this possibility, the couplings among 
the vector fields are chosen so that every quantum contribution disappears. In this way 
one is able to construct theories whose amplitudes are purely classical and, in principle, can 



generate topological terms describing interactions among an arbitrary number of random 
chains. Apart from a constant factor, each of these terms can be identified with the contri- 
bution of a tree-level Feynman diagram appearing in a non-Abelian C-S field theory ||26|| - 



p7| . However, it is exactly the freedom of choosing these factors that allows the fixing of the 
topological costraints. The above findings are used to build a model of entangling polymers 
in which the topological interactions induced by the Gauss linking invariant in the standard 
Edward approach are corrected by higher order topological interactions among three trajec- 
tories. The corresponding topological invariant has a simple physical interpretation which is 
expressed in terms of magnetic fields generated by fictitious charged particles moving along 
the trajectories of the polymers. A generalization to the case of four trajectories is also 
outlined. 

The material presented in this paper is divided as follows. In the next Section we explain 
the treatment of the N polymers problem of [§ based on the use of the Gauss linking number 
to distinguish the topological states. With respect to |^ some new investigations are made. 
For instance, the role of the C-S fields as propagators of the collective modes which are 



relevant in the topological entanglement |T^] is explored in details. In Section three that 



approach is extended to include also higher order topological interactions. To this purpose, 
abelian theories are defined that generate topological interactions among three and four 
loops. Finally, the Conclusions are drawn in section four. 

II. THE A^-POLYMERS PROBLEM 

Let Pi, ... , P/v be a set of topologically entangling random chains at thermal equilibrium. 
If the step length a of the segments composing the chains is very small, one can describe 
the chains as trajectories in the space parametrized by vectors rj(sj), i = 1,...,N and 
continuous parameters Si, . . . ,Sn such that: 

0<Si<Li ri(0) = r:,r,(L,)=r, (1) 



Li coincides with the contour length of Cj. In first approximation the topological constraints 
will be imposed exploiting the Gauss linking number: 

and requiring the conditions: 

xiCi,Cj) = rriij (3) 

where m^ = rriji, ma = are a set of topological numbers. 

To describe the statistical mechanics of the chains we define the configuration probability 
Cjmidr}, {L}] {r'}, 0). This function measures the probability that the trajectories Ci have 
extrema (in the open case) at the points r'- and r, or a fixed point (in the closed case) 
in r^ = Fj for i = 1, . . . , iV. Moreover, they should fulfill the topological conditions (|]). 
In our notations {m} denotes the n x n symmetric matrix of topological numbers, while 
{r} = Fi, . . . , Ttv, {L} = Li, . . . , Ltv etc. In the path integral approach one obtains the 
following expression of Gimjdr}, {L}; {r'}, 0) |^: 

G{m}{{r},{L};{r'},0)= I .../, exp {-(A + A.)} I] IT ^(x(C., C7,) - m,,) (4) 

j>i 

where 

A ^ I; g />?(..) (5) 

is the action of a free random walk and 

-4- = ^ E I d^^l ds'jv^/''\r,{s,) - r,(4)) (6) 

takes into account the excluded volume interactions. For convenience, we have put: 

'"ij = \ '"ij = '"ji (J) 

[ v%l2 for I ^ J 

where the v^j are coupling constants with the dimension of a volume. The Fourier trans- 
formed of Gimidr}, {L}; {r'}, 0) with respect to the parameters rriij is: 



r, fT 



■ N 



r', Jr'^ 



N-l N 



G{x}{{r}, {L}; {r'}, 0) = ^ ^ri(si) . . . j^ I)r^(sjv)exp {- (A + A. + Al)} (8) 
with 

^2L = ^EEx(a,c7,)A., (9) 

Following the approach of Edwards, one would like to transform the above path integral into 
a field theoretical problem. First of all, we treat the excluded volume interactions. To this 
purpose, we introduce A^ Gaussian scalar fields 0i(r), . . . , 0Ar(r), with action 

2 N 

^<t> = \T. jd'vH{v')-'rh (10) 

The fundamental identity which relates the excluded volume term to a field theory amplitude 
is: 



^-Aev 



j V(f)i . . . V(f)N exp <^ -A{^} - ^ E f dsi(f)i{ri{si)) \ (11) 



In the case of topological interactions one needs instead Chern-Simons fields A^*-, and B^l 
with action: 

N-l N 

Sc.. = -//rEi:A«.(VxB») (12) 

j>i 

The C-S theory will be quantized in the Landau gauge, where the fields are completely 
transverse. In the coupling with the random chains only the following linear combinations 
of fields are relevant: 

C'" = E^A»! (13) 

j=2 ^" 

A^ I. N~2 

C« = E J-^) + E A..Bg ^ = 2,...,N-l (14) 

j=3 ^^ i=i 



j>i j<i 



and 



N-l 
C(^) = E ^^NB% (15) 

1=1 



The topological term appearing in the configurational probability (P) can be rewritten as 
an amplitude of the above C-S field theory as follows: 

JVAVB exp|-z^cs-«E^'''C»(r(s.))rfr(s,)| =6--^^^ (16) 



where 

N 

u 



VA.VB = I n ^A«I)B[;.)) (17) 



FIG. 1. Graphical interpretations of the topological interactions between two random walks Cj 

(i) (i) 

and Cj mediated by the Chern-Simons fields A)-{,'B)A- 

Exploiting the identities (|1T]) and ([T6| ) the configurational probability (^ becomes: 

N 

G{A}({r}, {L}; {r'}, 0) = (11 ^(r,, L,; r^, O|0„ C«))^^^_^AmB} (18) 

1=1 

where 

G'(r„L,;r:,O|0i,C«) =|'"'Pr,(s,)e-/o''^''^^' (19) 

i 

and 

C, = ^r,2(s,) +20,(r,) -^r,(s,) • C«(r,(s,)) (20) 

As we see from eqs. (0) the trajectories Ci, . . . , Cat are completely decoupled before averag- 
ing over the auxiliary fields {(f)}, {A}, {B}. Formally, G(rj, L,; r^, O|0j, C*-*-*) is the evolution 
kernel of the random walk of a particle in an electromagnetic field (i0i,C'*^). Thus, it 
satisfies the pseudo-Schrodinger equation |T^ : 



d a 



^■+i(, 



G(r„L,;r:,O|0„C«) = 5(L,)5(r,-r:) 



pLi 6" 
The covariant derivatives D, appearing in the above equation are given by: 



(21) 



Di = V + iC'^^ 



1 = 1,. ..,N 



(22) 



It is now convenient to perform a Laplace transformation of (^(rj, Lj; r^, O|0i, C^*^) with 
respect to the length Li. 



G{vi,Y[-Zi\(t)i,C 



{i)\ 



dL,e-"»^'G'(r,,L,;r:,0|(/),,C 



{i)\ 



(23) 



Accordingly, we are now considering the Laplace transformed configurational probability: 

G{x}{{r}, {r'}, {z}) = J^ dL, . . . j^ dL^exp I - ^ z,L, | G{x}{{r}, {L}; {r'}, 0) (24) 

Since the order of the integrations over the auxiliary fields and the lengths Li of the trajec- 
tories can be permuted, we have: 

N 



G|,|({r},{r'},{4) = {l[Gi^^^^^■^^M^^C^'^))w,{AuB} 



(25) 



2=1 



The new variables Zi play the role of Boltzmann-like factors which govern the distribution 
lengths of the random chains. The advantage of having performed the Laplace transforma- 
tions is that G{ri, r^; -2i|0i, C*-*-') satisfies a stationary pseudo Schrodinger equation: 



[z,-H,]Giv„r',;z,\<Pi,C^'^) = 5ir,-T[) 



(26) 



Here the Hamiltonian Hi is given by: 



o 



(27) 



The solution of eq. ( p6D can now be expressed in terms of second quantized fields fields 

^*,iJi,t = l,...,N: 



G{r,, r[; z,\<P„ C«) = Zr' J V^,V^:4;,{r,)4j:{r[)e 



*('t-'^p--^M 



(28) 



where F[%lji\ represents the Ginzburg-Landau free energy of a superconductor in a fluctuating 
magnetic field: 

F[iJi]= jdh ^|D,V'«r+(^.+^0oiV'ir (29) 

and Zi is the partition function of the system: 

Z, = Jvij,Vi:*e-''^^'^ (30) 

The auxiliary fields (pi can be eliminated from the configurational probability (^) integrating 
them out. The integration over these fields is non-trivial due to the presence of the factors 
Z^^ in the second quantized expression of G'(rj, r^; 2;j|0j, C^*^) in (pH]), but can be made 
Gaussian by exploiting the identity pM]: 



Zr^ = lim Z^-i (31) 

In this way 

^(r^r^sz^^CW) == hm Zr-i /p^,P^;^,(r,)^;(r:)e-^['^'l (32) 

The above equation should be understood as follows: the right hand side is first computed 
supposing that the replica index Ui is an arbitrary positive integer and then one performs the 
analytic continuation of the result to the point rii = 0. Now G(rj, r^; Zj|0j, C'*^) is a product 
of rii path integrals. To each one we associate a set of replica fields ipl*, ai = 1, . . . , rij, so 
that it will be convenient to introduce the multiplets: 

^, = (^,^...,co (33) 

vl/* = (^f,...,^r) (34) 

and to rewrite G{ri, r^; Zi|0i, C*-*-*) as follows: 

G(r„r:;^,|0„C«) = lim f V%V^If*i;l{T,)i;*\r'^ e'^l^J (35) 



where 



F[<if, 



ai=2 



6 



U^t'l +{z^+^<PiM: 



a,i\2 



and 



dh 



\D^!A^ + {zi + i(l).-Mi 



(36) 






(37) 



According to the rephca method, one supposes that it is possible to commute the hmit of van- 
ishing rephca index with the path integrations over the auxihary fields. Thus, substituting 
eq. (^) in (pSD and performing the integration over the scalars 0j, which is now Gaussian, 
one obtains the final expression of the configurational probability G{A}({r}, {r'}, {z}): 

G{A}({r},{r'},{.}) = 



N 

hm / Vmvm* / PAPB n kKr.OV'f (r;.)] exp{-^Ga«..} 

.7 = 1 



(38) 



•Acauss is the free energy of the topologically linked random walk written in terms of C-S 
and second quantized fields. The subscript refers to the fact that the topological constraints 
have been imposed using the Gauss linking number. After a rescaling the complex scalar 
fields of the kind 



^i 



-^, 



\1/* 



-v]/* 



(39) 



where M is a mass parameter, the explicit expression of Acauss is: 



N 



Ac 



i=l 



Gauss = iScs + E / ^^r iDi^iP + m^l^ 



+ 



N 






/"t/^rl^iPl^jl 



(40) 



where 



m^i = 2Mzi (41) 

The above action describes a 0{n) model coupled to Chern-Simons fields in the limit n = 0. 
The analogous of the Planck constant is here the constant h = ^ and has been set equal 
to one in (|40|). The topological fields are not just auxiliary, but play a physical role, since 
they propagate the long-range interactions that impose the topological constraints (^. One 
may argue at this point that the number of C-S fields used in the present approach, which 
is A^(A^ "~ 1)) is highly redundant with respect to the physical number of degrees of freedom 
involved if A^ > 3. This number can be computed exploiting refs. [0, |2^, where a set of 



collective modes which are relevant in the topological interactions has been constructed in 
terms of the so-called bond vector densities: 



u,; r 



i dTiSir-r^) t = l,...,N (42) 



More precisely, the collective modes are linear combinations of the N bond vector densities 
in the Fourier space, where 

Ui(q) = / dviC^^-^^ i = l,...,N (43) 



Indeed, the Gauss linking invariant (0) can be expressed as follows: 

X(a, Q) = / ^ q • (u^(q) X u^(-q)) (44) 

J q 

As a consequence, if there is a number M < A^ of random chains which have non-trivial 
topological relations with the others, i. e. the maximum rank of the matrices {rriij} and 
{Xij} is M, there are at most M degrees of freedom to be propagated. Of course, this is 
not in contradiction with our result. As a matter of fact it is possible to see that, exploiting 
the equations of motions, the number of C-S fields in the action (^) can be reduced to A^. 
However, the C-S theory obtained in this way is not universal, since it cannot describe all the 
topological states of the system. The reason is that after the reduction the C-S propagators 
depend on the parameters X^j and become singular in the limit in which some of them vanish. 
In general, it has not been possible to build a suitable Abelian C-S field theory with less 

10 



than A^(A^ — 1) fields without encountering the problem of diverging propagators wherever 
rank[X\ < N or without resorting to a complicated parameterization of the matrix {A} 
provided for instance by the solution of the following algebraic system of equations: 

N 

^ij = Yl VikVkj (45) 

fc=l 

III. INCLUDING HIGHER LOOP INTERACTIONS 

The Gauss linking number describes a topological interaction between two loops and it 
is quite a poor topological invariants. Thus it would be interesting to include in the above 
approach also higher order topological interactions. To begin, we consider an interaction 
r3(Cj, Cj, Ck) among three loops, where i < j < k. To determine T^i^Ci, Cj, Ck) we construct 
a suitable topological field theory with action: 

53(^, J, k) = e^^^ J d'xa%,^d^a%,^ + e^^^ J d^ xh%,^djb%,^ 

^fiup J d ^(^^ji(ij)9uCp(,ij-^ + K^^e^yp j d a^a^(jfc)^iy(ifc)'^p(ii) (46) 

where the fields %(jfc)' • • • ' ^Lfc) ^^^ purely transverse. Ss{i,j,k) describes at the classical 
level a topological field theory which is not gauge invariant. Moreover, it is easy to convince 
oneself that the theory has no radiative corrections that could spoil its topological properties. 
Despite of this fact, there are nontrivial amplitudes as for instance the following correlator: 

G^. (a, C„ Ck) = {e'^c. ^^^.^^^^'^c, <^.,^^^)n^ §,^ <^l^,,(r.)^ ^^^^ 

The above amplitude can be exactly computed and one obtains: 

G^.ja, C,, Ck) = exp {A},r(a, C,, Ck)] (48) 



where 



T{C„C„Ck)=i drfi rfrf/ rfr^/,^^(r„ r„ r^) (49) 

JCi J Cj J Ck 
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and 

/«^^(ri,rj,rfe) = e'^'P J dhG ^air - ri)G^(i{r - rj)Gp.,{r - rt) (50) 

In the above equation 

'^ '^ |ri -r2|^ 

r'^ being the components of the vector r. An analogous of r(Cj, Cj, Ck) in the case C, = Cj = 
Ck has been studied in connection with perturbative calculations of self-linking invariants 



in non-Abelian C-S field theories. After a volume integration in (pOD, the right hand side of 



eq. (pop has a complicated expression, which has been evaluated in |^. With respect to 
P7| , however, one does not need path ordering of the trajectories, so that it is possible to 
perform the volume integration using a different strategy. To this purpose, let us define the 
currents: 

j^,){r)^f drf6{r-ri) (52) 



and the magnetic fields: 



i3(o,(r) = 6,., / dr'^T—^ (53) 



with V-i3(i) = and V x B(^i) = 47rj(;-). The vector potentials corresponding to these magnetic 
fields are: 

A,,(r)=i ^^fa^r'hl^. (54) 

Jci \r — r/l'^ J \r — r | 

ad satisfy the relations 

V ■ ^(0 = V X A^i) = B^i) (55) 



A^(,) = -47rj(,) (56) 

In terms of the magnetic fields B^i-. we have: 

12 



r(a,C,,C,) = jdhB^c^iv) ■ 4)(r) X i3f,^(r) (57) 

The space integral can be eliminated using the Stokes theorem and after some calculations 
one obtains: 

-^r(C„ C„ Ck) = - f f dSk ■ %(rfe) X 4)(rfc) + / / dS, ■ %)(r,) x %)(r,) 
47r J •iT.Cf, -I •''^Cj 



rfS,-%)(r,,)x%)(r,) (58) 

where the Sc-., Sc^-, Sc^, are arbitrary surfaces having respectively Cj, Cj, C^ as borders and 
infinitesimal surface elements dSi,dSj,dSk- To include the interactions among the loop 
trajectories, we extend the previous configurational probability of eq. (§) as follows: 

G{m,M}i{^}AL};{r'},0)= I ...I exp{-(A + A.)} 

N-1 N N-2N-1 N 

n n ^(xiQ, c,) - m,,) n n n '^(^(a, c,, o - m,,,) (59) 

i=l J=2 i=l J = 2 fc=3 

j>i j>i k>j 

After performing the Fourier transformations with respect to the topological numbers mij 
and Mijk one obtains: 

G|AA}({r}, {L}; {r'}, 0) = / ^ 2^ri(si) . . . / ^ Prjv(siv)exp {- (A + A. + Al + Al)} 

(60) 

with 

AT 

A,. = z Y. {A}.r(a,C7,,Cfc)} (61) 



i<j<k 



At this point we rewrite the higher loop topological interactions appearing in ( p(]| ) as a field 
theory amphtude using eqs. (^) and (^Of ): 



r'[a]P[6]I?[c]r'[S]r'[6]r'[c]e-^^3^ 
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fj e'^^O, ^r?a«„(r.)^^^/e, <^^,l.,(r.)^n3 /^^ ^r^c«^,(r,) 



i.j'jfc — 1 
i<j<k 



o-iAzL 



(62) 



.» 



;;(fc) 



X>[a], . . . , "Dfc] denotes the measure over the C-S fields O'Jjk)^ • • • > c|,'(^j) and 

N 



(63) 



i.j.fc — 1 
i<j<k 



Again, the C-S fields decouple the topological interactions. The integration over the tra- 
jectories of the random chains can now be performed following the same strategy of the 
previous section. In this way one obtains the following expression of the configurational 
probability in terms of the Laplace variables Zi, . . . , z^: 

G{A,A}(W,{r},M) = 



where 



lim Vmvm* / VAVBV\a] . . . V\c] FT \i)l{Yi)i)f{ri)' 

ni,...,n]v^O J 1 '- 

1=1 



^-Ah 



N 



Ah = iScs + ^S^L + E / ^'r [|Df il/.p + ml\m 

8=1'' 



+ 



(64) 



E— r^y^'r|v^,n^,r 

The covariant derivatives that include the three loop interactions are given by: 



D 



--*^ V + i(C(*) + d^'^) 



(65) 



(66) 



with 



N N N 



s<j<k 



i,k^l 
i<s<fc 



i<j<s 



m 



(67) 



We notice that within the above approach it is possible to include also topological interac- 
tions among four or more trajectories. Let us consider for instance a topological four loop 
interaction T{Ci, C2, C3, C4). To generate such interaction we define the following action: 
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^4(1,2,3,4) = ^^|rf3^eM^p^»5,6W + 



Ai J d^xe^-'^Pa^^^'^afaf + A2 J d'xt''^"afhfaf (68) 

As in the previous case, the theory is topological but has no gauge invariance. This can be 
dangerous, since radiative corrections may arise which spoil the topological properties of the 
theory, but it is easy to see that the above theory has no quantum contributions. 
The relevant correlation function to be considered here is: 



GA,A,(7i,72,73,74) = (ne-^-/<^^l^^) 



(69) 



i=l 



The above amplitude can be exactly computed and the result is: 



G'aiA2(7i, 72, 73, 74) = exp <^ -iAiAs 



n7. 

i=l 



r(Ci, C2, C3, C4) 



(70) 



where 



r((:7i,C2,C3,C4) = n [/ dxt'] I <fx f d'^ye^^'e^''^ 



X Gy,„{x - y)Gxai{x - Xi)Gua2{x - X2)Gpa^{y - X-i)Grai{y - X4) 



(71) 



Eq. ( [TTD describes a topological interaction among four loops. Unfortunately, the elimination 
of the double volume integral has not been possible following the strategy of Section II. For 
this reason, an expression of r(Ci, C2, C3, C4) in terms of the magnetic fields (^) could not 
be derived. 



IV. CONCLUSIONS 

In the first part of this work the statistical mechanical problem of a system of A^ polymers 
whose topological interactions are governed by the Gauss linking number has been mapped 
to a field theory following ref. [^. In the model obtained in this way the C-S fields play a 
physical role, since they mediate the topological forces which impose the constraints (|^). We 
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notice that the relation ([I6D, which expresses the topological contributions appearing in the 
first quantized version of the configurational probability (^) in terms of field amplitudes, can 
be reproduced also by means of other Abelian C-S field theories. Each of these theories differs 
from the other by the number of fields. However, the requirement that their propagators 
should not diverge when rank[X\ < N implies that the independent fields should be at least 
A^(A^ — 1). On the other side, if one starts with N > N{N — 1) C-S fields, it is always 
possible to reduce their number to A^(A^ — 1) by exploiting the equations of motion. For this 



reason, the model of topologically entangling polymers given by eqs. (p8|)-(p]) is unique. 

In Section II the field theories are essential in order to decouple the random chains 
and to rewrite the configurational probability in a second quantized form. However, they 
start to play an even more active role in Section HI. In fact, here the field amplitudes are 
fundamental to provide the explicit expression of the topological interactions among three or 
more polymer trajectories. In the case of a three loop interaction, the topological invariant 
r{Ci,Cj,Ck) has also a nice physical interpretation in terms of magnetic fields given by 
eq. (^) as the Gauss linking invariant. With respect to non- Abelian C-S field theories, 
the advantage of the theories of transverse vector fields defined here is the possibility of 
generating single topological invariants r{Ci,Cj,Ck) without the problem of spurious non- 
topological contributions or of higher order corrections. Moreover, the freedom to choose 
the parameters A^^/^ is crucial in order to impose constraints on the T{Ci,Cj,Ck) as shown 
by eqs. (H) and {^. 

Concluding, the field theoretical approach illustrated in this paper solves in principle the 
most serious drawback present in the analytical approach of Edwards, i. e. the use of the 
Gauss linking invariant to specify the topological status of a system of random chains. How- 
ever, more work and, possibly, numerical simulations are still needed in order to evaluate 
the phenomenological implications of the new topological terms introduced in the configu- 
rational probability (|64D and to make contact with experiments. 
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